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Bilayer silicene has richer physical properties than bilayer graphene due to its buckled structure together with 
its trigonal symmetric structure. The trigonal symmetry originates in a particular way of hopping between two 
silicenes. It is a topologically trivial insulator since it carries a trivial Z2 topological charge. Nevertheless, its 
physical properties are more akin to those of a topological insulator than those of a band insulator. Indeed, 
a bilayer silicene nanoribbon has edge modes which are almost gapless and helical. We may call it a quasi- 
topological insulator. An important observation is that the band structure is controllable by applying the electric 
field to a bilayer silicene sheet. We investigate the energy spectrum of bilayer silicene under electric field. Just 
as monolayer silicene undergoes a phase transition from a topological insulator to a band insulator at a certain 
electric field, bilayer silicene makes a transition from a quasi-topological insulator to a band insulator beyond 
a certain critical field. Bilayer silicene is a metal while monolayer silicene is a semimetal at the critical field. 
Furthermore we find that there are several critical electric fields where the gap closes due to the trigonal warping 
effect in bilayer silicene. 
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INTRODUCTION 

Silicene is a monolayer of silicon atoms forming a two- 
dimensional honeycomb lattice UHl]. It has a relatively large 
spin-orbit (SO) gap, and its intrinsic property is the buckled 
structure ||4|-|6|] owing to a large ionic radius of silicon. Sil- 
icene has richer physical properties than graphene due to this 
property. First of all, it is a topological insulator, character- 
ized by a full insulating gap in the bulk and helical gapless 
edges 1 6]. Furthermore, the band structure is controllable by 
applying the electric field E z to a silicene sheetj7|]: It has 
been shown that silicene undergoes a topological phase tran- 
sition from a topological insulator to a band insulator as \E Z \ 
increases. It is a semimetal at the critical field E CY due to its 
linear dispersion relation. It has many attractive and remark- 
able properties [fsl-fioh . Silicene may be the most promising 
material now available as a topological insulator flllUl4ll . 

In this paper we analyze the band structure of bilayer sil- 
icene, which was manufactured! 15] experimentally only re- 
cently. Though bilayer silicene has a trivial Z2 topological 
charge, we show that its physical properties are more akin to 
those of a topological insulator than those of a band insula- 
tor. When we switch off the Rashba SO interaction (Ar = 0) 
and the interlayer SO interaction (Ai nte r = 0), bilayer silicene 
shares the same edge-mode properties with a topological in- 
sulator: It has a full insulating gap in the bulk and helical gap- 
less edges. However, helical gapless edge modes are not topo- 
logically protected because of a possible mixing between the 
two gapless modes present in each edge. Indeed, a small gap 
opens when these interactions are taken into account. Never- 
theless, the physical properties are very much similar to those 
of a topological insulator. This is particularly so under electric 
field, as we shall soon see. It would be reasonable to call such 
a system a quasi-topological insulator. 

Here we recall a recent proposal of bilayer graphene [fl6l 
with rather unphysical parameters for study of topological in- 
sulator. The model Hamiltonian is very similar to ours, and 
there are much in common. Nevertheless, the great merit of 



our theory is that silicene is a realistic material, and any ex- 
perimental test on our results is feasible. Furthermore, it is 
remarkable that the band structure is controllable by applying 
the electric field to bilayer silicene. Just as monolayer silicene 
undergoes a phase transition from a topological insulator to 
a band insulator at a certain electric field[7], bilayer silicene 
makes a transition from a quasi-topological insulator to a sim- 
ple band insulator beyond a certain critical field. It is a metal 
at the critical field. There exists also a new feature. We find 
that the band gap closes a few times when we increase electric 
field, which is absent in the case of monolayer silicene. It is a 
trigonal warping effect. 

Bilayer silicene has an additional degree of freedom how to 
stack two buckled monolayer silicenes. There are four types 
of stacking even if we concentrate on the Bernal (AB) stacking 
(Fig.l). For example, the forward stacking has an electron- 
hole symmetry but the backward stacking does not. It is pos- 
sible to determine the type of stacking of a sample by experi- 
mental measurement. 

This paper is organized as follows. In Section II we intro- 
duce the Hamiltonian for bilayer silicene and study the band 
structure of a nanoribbon. We point out by numerical calcula- 
tion that the edge modes are almost gapless and helical though 
bilayer silicene has a trivial Z2 topological charge. In Section 
III we derive the effective Dirac theory to describe the low- 
energy physics around the K and K' points. In Section IV 
we investigate the band structure under homogeneous electric 
field based on the effective Dirac theory. Bilayer silicene is 
shown to become metallic at the critical electric field, though 
the band structure depends considerably on the type of stack- 
ing of two buckled silicenes. 



SILICENE AND BILAYER SILICENE 

Silicene consists of a honeycomb lattice of silicon atoms 
with two sublattices made of A sites and B sites. Due to the 
buckled structure the two sublattice planes are separated by a 
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FIG. 1 : (Color online) Schematic illustration of AB stacking bilayer 
honeycomb lattice. Bond connecting sites A, B in the top layer is 
indicated by solid line, while bond connecting sites A, B in the bot- 
tom layer by dashed line. The interlalyer interactions are indicated 
by arrows together with the couplings t±, £3 and A^ter- 



distance, which we denote by 21 with i — 0.23 A. The states 
near the Fermi energy are 7r orbitals residing near the K and 
K' points at opposite corners of the hexagonal Brillouin zone. 
We refer to the K or K' point also as the point with the 
valley index 77 = ±1. 

The monolayer silicene system is described by the four- 
band second-nearest-neighbor tight binding model [|17ll . 
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where c\ a creates an electron with spin polarization a at site 
i, and (i, j) / run over all the nearest/second-nearest 

neighbor hopping sites. The first term represents the usual 
nearest-neighbor hopping with the transfer energy t — 1.6eV. 
The second term represents the effective SO coupling with 
Aso = 3.9meV, where a = (a x , (j yi a z ) is the Pauli matrix of 
spin, and = +1 if the second-nearest-neighboring hopping 
is anticlockwise and vij = — 1 if it is clockwise with respect 
to the positive z axis. The third term represents the Rashba 
SO coupling with Ar = 0.7meV, where fii = ±1 for the A 
(B) site, and d^ = d^j \d^ \ with d^ the vector connecting 
two sites i and j in the same sublattice. Monolayer silicene 
has been shown to be a topological insulator|0|7|]. 

We model a bilayer silicene as two coupled buckled hexag- 
onal lattices including inequivalent sites A, B in the top layer 
and A, B in the bottom layer, respectively. These are arranged 
according to Bernal (A- B) stacking, as illustrated in FigQ] 
There are four types of Bernal stacking in bilayer silicene, 
since the top and bottom layers can be buckled upward and 
downward independently. Irrespective of the type, the bilayer 
silicene system is described by the eight-band second-nearest- 
neighbor tight binding model, 
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where H$^ is the Hamiltonian (Q]) for the top (T) or bot- 
tom (B) layer; i^inter is the interlay er Hamiltonian, where the 
first term is the nearest-neighbor interlay er vertical hopping 
between A sites of the top layer and B sites of the bottom 
layer with the coupling t±, and the second term is the next- 
nearest-neighbor interlayer hopping between B sites of the 
top layer and A sites of the bottom layer with the coupling £3. 
There are three ts hopping vectors from one site, as implies 
the trigonal symmetry of the bilayer system. The parameters 
t± and ts depend on the type of Bernal stacking, but we ex- 
pect t± > ts ~ O.lt in all of them. The third term describes 
the interlayer SO interaction [fl8l [l9ll with the coupling Afr- 
its magnitude would be the same order as that in graphene 
(Ainter ~ 0.5meV). 

It is readily shown that the band structure of the Hamilto- 
nian d21) is gapped, and hence bilayer silicene is an insulator. 
To study what type of insulator it is, we have analyzed the 
energy bands of a bilayer silicene nanoribbon based on the 
tight-binding model (EJ), whose result we give in Figj2j The 
band structure consists of the bulk modes and the edge modes, 
where the bulk modes are gapped. There exists an intriguing 
feature with respect to the edge modes. 

When we switch off the Rashba and interlayer interactions 
(Ar = A^ter = 0) there emerge eight helical zero-energy 
modes in the edges [FigEJa)]. This is what we expect since 
the number is simply twice as much as that in monolayer 
silicene |@, 0]. Namely, four helical zero-energy modes arise 
from each monolayer silicene. Indeed, according to the bulk- 
boundary correspondence rule, and since the Chern number of 
the bilayer system is shown to be twice that of the monolayer, 
the number of helical edge states should also be double. 

However, when we switch on these interactions (ArAi^ 7^ 
0), a spin mixing occurs, as turns a two-band crossing into 
a two-band anticrossing and opens a small gap for the edge 
modes [FigEJb)]. This makes a sharp contrast to the mono- 
layer case, where the gapless edge modes are protected topo- 
logically against perturbations. Despite each layer behaves as 
a topological insulator, an even number of such layers renders 
the overall system topologically trivial that is characterized by 
a vanishing Z2 charge. The topological triviality of a bilayer 
system may be understood as follows. In the bilayer system, 
backscattering between channels at the same edge moving in 
opposite directions with opposite spins is not forbidden, in 
contrast to the case for the monolayer system. We note that a 
similar behavior has been pointed out in bilayer graphene [ll6ll . 

Nevertheless, the properties of bilayer silicene are more 
akin to those of a topological insulator than those of a band 
insulator. When Ar = Ai nte r = the edge modes are purely 
helical, which is one of the characteristic feature of a topo- 
logical insulator. When ARAi nte r / a spin mixing occurs 
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FIG. 2: (Color online) One-dimensional band structure for a bilayer 
silicene nanoribbon. It consists of the bulk modes and the edge 
modes. The vertical axis is the energy gap A in unit of t. The hor- 
izontal axis is k x . (a) When the Rashba and interlayer interactions 
are switched off (Ar = Ai nte r = 0), there are eight bands crossing 
the gap, which are edge modes. There are four due to the left (L) 
and right (R) edges of the top (T) and bottom (B) silicenes, with 
each of them doubly degenerated with respect to the spin t and > | r . 
The gapless edge modes are helical, (b) When these interactions are 
switched on (ARAinter / 0), a crossing turns out to be an anticrossing 
due to a spin mixing, as induces a small gap. The edge modes are 
almost helical except for the anticrossing point. In fact, spin mixing 
occurs only in the vicinity of the anticrossing point. We have taken 
Ar/Aso = 1/2 and Ainter/ Aso = 1/4 for illustration to emphasize 
the gap. 



only around the anticrossing point, and the edge modes are 
almost helical away from the point, as we can confirm numer- 
ically. As we shall soon show, by applying a critical electric 
field, the gap becomes precisely zero at the K and K' points: 
Bilayer silicene becomes a metal due to a parabolic disper- 
sion, just as monolayer silicene becomes a semimetalj7|] due 
to a linear dispersion. It is reasonable to call such an object a 
quasi-topological insulator. 



LOW-ENERGY DIRAC THEORY 

We analyze the physics of electrons near the Fermi energy 
based on the low-energy effective Hamiltonian derived from 
the tight binding model ©. It is described by the Dirac theory 
around the point as 
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FIG. 3: (Color online) The band structure of bilayer silicene near 
the K point based on the Dirac Hamiltonian (4j). The vertical axis 
is the energy gap A in unit of t. The horizontal axis is k x . It is 
clearly gapped. There is no symmetry for k x < — > —k x due to the £3 
hopping term, reflecting the trigonal warping. 
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where r^ B (ri ayer ) is the Pauli matrix of the sublattice (layer) 
pseudospin, = ^at = 5.5 x 10 5 m/s is the Fermi velocity, 
and a = 3. 86 A is the lattice constant. We show the band 
structure in Figj3] 

The Hamiltonian contains many parameters. Actually, Ar 
and Ai nter are small constants with respect to the other param- 
eters. Hence, it is a good approximation to set Ar = Ai nter = 
in most cases. In what follow, although we carry out numeri- 
cal analysis by including the nonzero effects of Ar and Ai nte r, 
we develop an analytic formalism by assuming Ar = Ai nte r = 
0. It enables us to make a simple and clear physical picture on 
the basis of analytic formulas. For instance, the spin s z = ±1 
becomes a good quantum number in this simplification, where 
the spin-chirality is given by s = rjs z . We can always check 
the validity of the approximation numerically. 

We rewrite down the Hamiltonian H+ explicitly in the basis 
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while the off-diagonal blocks are 



Case 1 



Case 2 



H12 = H 2 i = 



hv?k- 



(11) 



hv?k+ 

with vs = ^ats and k± = k x ± ik y . The diagonal elements 



are 



(12) 



where ti ayer = ±1 for the upper (lower) silicene sheet. We 
have introduced this parameter for a later convenience, though 
it is irrelevant here. 

We make a further simplification of the theory by reducing 
the 4x4 Hamiltonian to the 2x2 Hamiltonian, which is valid 
as the effective two-band Hamiltonian for the lower energy 
bands[20]. Provided ts < t±, the reduced Hamiltonian is 
given by 



where G 2 2 is me 2x2 Green function, 
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with the energy e from the Fermi level. The effective Hamil- 
tonian around the point reads 



H\ =h B + h w - sXso 
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where we have used \e\ <C t±. The Hamiltonian Kb describes 
the quadratic Dirac Hamiltonian, which is familiar in bilayer 
graphene. The Hamiltonian h w describes the trigonal warping 
effects. We note that the reduction into the 2x2 Hamiltonian 
is not available when t 3 > t± . 



BILAYER SILICENE IN ELECTRIC FIELD 

We take a bilayer silicene sheet on the xy -plane, and ap- 
ply the electric field E z perpendicular to the plane. It gener- 
ates a staggered sublattice potential oc ±21E Z between silicon 
atoms at A sites and B sites in a single silicene sheet, and also 
a staggered sublattice potential oc ±2LE Z between two sil- 
icene sheets. How they enter into the Hamiltonian depends on 
the type of four Bernal stacking. They are 

1) The top layer buckles downwardly and the bottom layer 
buckles downwardly (Fig|4)). 
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FIG. 4: (Color online) Illustration of the forward Bernal (AB) stack- 
ing bilayer silicene. The order of sites is B,A,B,A for the case 1 
and A, B, A, B for the case 2 from the top to the bottom. 



2) The top layer buckles upwardly and the bottom layer 
buckles upwardly (Fig@]). 

3) The top layer buckles upwardly and the bottom layer 
buckles downwardly (Fig 2]). 

4) The top layer buckles downwardly and the bottom layer 
buckles upwardly (Figj5j). 

Here, upward (downward) means that A site is higher (lower) 
than B on the same layer. 

Silicon atoms prefer having the neighboring atoms in reg- 
ular tetrahedron directions, and hence the case 1 seems most 
likely. The cases 3 and 4 are identical by reversing the whole 
system, and the energy spectrum is identical between them. 
Thus it is enough to investigate the cases 1, 2 and 3. These 
three systems can not be transformed only by the transla- 
tion. It is necessary to rotate 120 degrees after translation. 
We call the case 1 and 2 (3 and 4) as the forward (backward) 
Bernal stacking. The forward (backward) Bernal stacking bi- 
layer silicene is shown in Fig]?]©, where the order of sites is 
A,B,A,B (B,A,A,B). The electric energy depends on the 
order. The magnitude of the interlay er SO interaction Ai nte r 
depends on the stacking type. 

The effective Hamiltonian for the forward stacking is given 

by 

H = H BL + {Lr ] r r + tr™) £ ViK<L^ (18) 

ia, 

while the one for the backward stacking is given by 

H = H BL + (L- It™)t^ E ViKclaa, (19) 

ia. 

where Hbl is given by ©. Accordingly, the Dirac Hamilto- 
nian dU) is modified as 

H v = H° v +E(s,t™,t l D, (20) 
where H® is the Dirac Hamiltonian ©, and 

E(s, if , 4? er ) = -sAsoif + (£r^ yer + tr^)E z (21) 
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FIG. 5: (Color online) Illustration of the backward Bernal (AB) 
stacking bilayer silicene. The order of sites is B,A,A,B for the 
case 3 and A, B, B, A for the case 4 from the top to the bottom. 



for the forward stacking, and 
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for the backward stacking. The modification is very simple for 
the four-band theory: The electric field is incorporated into the 
Dirac Hamiltonian ([8]) only by changing the diagonal elements 
E(s, d ayer ) in ® and (Ql) with (ED or d22]). 

We first investigate the case 1, and then the cases 3 and 
4. The case 2 is investigated at the end since the two-band 
Hamiltonian is unavailable. 



Forward Bernal Stacking in Electric Field: Case 1, BA-BA type 
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FIG. 6: (Color online) (a) Band gap A of the forward stacking sil- 
icene as a function of the electric field E z for the case 1. The 
gap closes at the critical points E CI o, E CI i, E ci2 , E cr 3. Bilayer sil- 
icene is an insulator except for these points. We note that the crit- 
ical points E cr2 and E cr s are missed in the effective two-band the- 
ory, (b) The band gap in the limit £3 — > 0, where all the critical 
points collap to one. Solid (dashed) curve is the band gap obtained 
based on the Dirac Hamiltonian d20t by taking Ar/Aso = 1/2, 
Ainter/Aso = 1/4 and Xso/t = 1/10 for illustration to emphasize 
the effects (A = Ainter = for comparison). They show qualitatively 
the same behavior. 



In the case 1 of the forward Bernal stacking, the effective 
two-band Hamiltonian under homogeneous electric field is de- 
rived as 

H eff =h B + h w + ra D ^ q _^ ^ 

(23) 

with the Dirac mass 

m D = (L - £) E z - s\ so . (24) 

The energy spectrum has the electron-hole symmetry. The 
spectrum becomes gapless when mo = 0, or E z = E cy q with 

which occurs at k± =0. We can show that the band gap is 
given as 

A = 2|m D |, (26) 



for IE 1 1 < E CY o. See also Fig J6] which we have derived by nu- 
merical calculations, where we compare the two results with 
Ar A^ter ^ and with Ar = Amter = 0. We show the band 
structure at the critical electric field in FigJ7](a). The band 
structure at the Fermi energy is parabolic, which is reminis- 
cence of bilayer graphene. We conclude that the system is 
a quasi-topological (band) insulator when the Dirac mass is 
negative (positive). 

The gap closes at E z = ±E CY o, where it is a metal due 
to gapless modes exhibiting a parabolic dispersion relation: 
See the band structure at E z = E CY o in FiglTJa). It follows 
that up-spin (s z = +1) electrons are gapless at the K point 
(rj = +1), while down-spin (s z = — 1) electrons are gapless at 
the K' point (77 = — 1). Namely, spins are perfectly up (down) 
polarized at the K (K') point under the uniform electric field 
E z = E cy q. Recall that monolayer silicene is a semimetal due 
to a linear dispersion at the critical point[|7|] . 

There exists one more critical point implied by the 2 x 2 
Hamiltonian, where the spectrum become gapless. The con- 
dition is that the Hamiltonian (f23l) itself becomes zero, as im- 
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FIG. 7: (Color online) Band structure of the forward Bernal (AB) 
stacking bilayer silicene (a) at the critical field E cr o and (b) at the 
critical field E CT i for the case 1 . The horizontal axis is k x . (a) The 
gapless mode is parabolic, showing that it is a metal at E = E cr o. 
(b) At E — .Ecri, there is no symmetry for k x < — > —k x due to the 
£3 hopping term, reflecting the trigonal warping. A 3-dimensional 
figure is given in Fig(8] 




FIG. 8: (Color online) Band structure of the forward Bernal (AB) 
stacking bilayer silicene at the critical field E cr i for the case 1. The 
Dirac cone is split into three cones, where the band touches the Fermi 
level. We have taken Ar/Aso = 1/2, Ainter/Aso = 1/4 and Aso/£ = 
1/10 for illustration. 



plies 



h 2 vl 



(L - t)E z - sXso ~ {(£ + I) E z + sAso} k_k + = 0, 
t ± 

hv 3 k+ - '-^±k 2 _ = 0. 



t± 



(27) 
(28) 



There are two solutions with k y = 0. The first one is k x = 
0, where E z is given by d25t . The new solution is given by 
solving 
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as illustrated in Fig|3b). There are two more solutions by 
performing ±2tt/ 3 rotations in the (k x ,k y ) plane around (0, 0) 
with the same critical field E CY i due to the trigonal symmetry, 
as illustrated in Fig HI 

Actually there is two more critical point, E CY 2 and E CT s, 
which are present in the original Hamiltonian (|8]> with (1211 
but missed in the reduced two-band model (l23k as we demon- 
strate in FigJSJa). Note that E cri — » E cy q for i = 1, 2, 3 as 
£3—^0. Namely, these three critical points are generated by 
the £3 hopping effect. 



Backward Bernal Stacking in Electric Field: Cases 3 and 4, 
BA-AB and AB-BA types 

The cases 3 and 4 are identical by reversing the whole sys- 
tem. In the case of the backward Bernal stacking, the effec- 
tive two-band Hamiltonian under homogeneous electric field 
is derived as 



i/IL =h B + h w +m D 



1 
-1 



£E Z 



1 
1 



(L^ + .Aso) / _ k k _ 



&4 
t 2 



-£E Z 



k-k + 
k+k- 



with the Dirac mass 



ra D = LE Z - s\ so . 



(31) 



(32) 



The energy spectrum has no electron-hole symmetry. The 
Fermi energy is moved to —IE Z by the presence of the electric 
field. The spectrum becomes gapless at E z = E cy q with 
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which occurs at k± = 0. We can show that the band gap is 
given as 



A = 2X S0 -2L\E Z \ 



(34) 
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FIG. 9: (Color online) (a) The band gap A of the backward stacking 
silicene as a function of the electric field E z for the cases 3 and 4. 
The gap closes at the critical points E cr o, E a i, E cr 2 and E CV 3. Bi- 
layer silicene is an insulator except for these points. We note that 
the critical points E cr 2 and E cr 3 are missed in the effective two-band 
theory, (b) The band gap in the limit £3 — ► 0, where all the criti- 
cal points collap to one. Solid (dashed) curve is the band gap ob- 
tained based on the Dirac Hamiltonian (f2(Tb by taking Ar / Aso = 1/2, 
Ainter/ Aso = 1/4 and Aso /t = 1/10 for illustration to emphasize the 
effects (A = Ainter = for comparison). The band gap is identical 
for the cases 3 and 4. 



for \E\ < E cy q. See also Fig(9]we derive by numerical calcula- 
tions.. We note that the energy spectrum is identical between 
the cases 3 and 4. We show the band structure at the critical 
electric field in FigfTUfa). The band structure at the Fermi en- 
ergy is parabolic, which is reminiscence of bilayer graphene. 

There exists one more critical field implied by the 2x2 
Hamiltonian, where the spectrum become gapless. The con- 
dition is that the Hamiltonian (l23l) itself becomes zero, as im- 
plies 



h 2 v 2 

LE Z - sXso - — 2^ (LE Z + sA S0 ) k_k+ = 0, 
1 1 
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FIG. 10: (Color online) Band structure of the backward Bernal (AB) 
stacking bilayer silicene at the critical fields E CT o and E cr i for the 
cases (3) and (4). The gap structure A for the case 4 is given by — A 
for the case 3. The gap closes at k x = for E cr o, but at k x ^ 
for E cr i due to the trigonal warping. The Fermi level is the one at 
E z = 0. 
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There are two more solutions by performing ±27r/3 rotations 
in the (k x ,k y ) plane around (0, 0) with the same critical field 
E CY i due to the trigonal symmetry. 

Actually there are two more critical points, E cv2 and E cr s, 
which are present in the original Hamiltonian © with (l22t 
but missed in the reduced two-band model (I3TK as we demon- 
strate in Fig|9] It is intriguing that the band is flat for E z > 
E CV 2. Note that E CY i — >> E CT o but E cr2 — » oo as ts — >> 0. 
Namely, these two critical points are generated by the ts hop- 
ping effect. 



Forward Bernal Stacking in Electric Field: Case 2, AB-AB type 

Finally we investigate the case 2. In this case, since we 
have ts > t±, the two-band Hamiltonian of the type (IT3l) is 
unavailable. We are unable to present analytic formulas. We 
have calculated the band gap numerically, which we show in 
FigHUThe number of the critical electric field is found to be 
reduced to two. 



There are two solutions with k v 



0. The first one is k r 



0, where E z is given by d33t . The new solution is given by 
solving 



± 



t ± y/LE z - sXso 
hv?\/LE z + s Aso 



V3t± 



(37) 



DISCUSSIONS 

Bilayer silicene has such a peculiar feature that its physi- 
cal properties are more akin to those of a topological insulator 
than those of a band insulator though it has a trivial Z2 topo- 
logical charge. Namely it is characterized by a full insulating 
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FIG. 1 1 : (Color online) (a) The band gap A of the forward stacking 
bilayer silicene as a function of the electric field E z for the case 2. 
The gap closes at the critical points E cx q and E CT i . Bilayer silicene is 
an insulator except for these points. 

gap in the bulk and almost helical gapless edges. It exhibits a 
similar behavior as monolayer silicene under electric field E z . 
Physical properties are very similar to those of a topological 
insulator for \E Z \ < E cy q, while they become those of a band 
insulator for \E Z \ > E cy q with a transition point E cy q. The gap 
closes precisely at the transition point, where bilayer silicene 
is a metal while monolayer silicene is a semimetal. 

As we have discussed in the case of monolayer silicene 0], 
therefore, by applying an inhomogeneous electric field, we 
would be able to create metallic regions anywhere within a bi- 
layer silicene sheet where helical zero-energy modes are con- 
fined. Furthermore, by rolling up a sheet, we may consider 
a double-wall silicon-nanotube, and by applying an electric 
field perpendicular to the tube axis, we would be able to create 
metallic channels parallel to the tube axis where spin currents 
are conveyed lllOll . 

I am very much grateful to N. Nagaosa for many fruitful 
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by Grants-in-Aid for Scientific Research from the Ministry of 
Education, Science, Sports and Culture No. 22740196. 
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